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The Casimir effect is a general phenomenon in physics, which arises when the vacuum fluctuation 
of an arbitrary field is modified by static or slowly varying boundary. However, its spin version is 
rarely addressed, mainly due to the fact that a macroscopic boundary in quantum spin systems is 
hard to define. In this article, we explore the spin Casimir effect induced by the zero-point fluctuation 
of spin waves in a general non-collinear ordered quantum antiferromagnet. This spin Casimir effect 
results in a spin torque between local spins and further causes various singular and divergent results 
in the framework of spin-wave theory, which invalidate the standard I/S' expansion procedure. Based 
on the spin Casimir torque interpretation, we develop a spin-wave expansion approach named as 
torque equilibrium spin wave theory (TESWT). In this approach, the spin Casimir effect is treated in 
a self-consistent way, and the spin-wave expansion results are free from singularities and divergences. 
A detailed spin-wave analysis of the antiferromagnetic spin-1/2 Heisenberg model on a spatially 
anisotropic triangular lattice is undertaken within our approach. Our results indicate that the 
spiral order is only stable in the region 0.5 < a < 1.2, where a is the ratio of the coupling constants. 
In addition, the instability in the region 1.2 < a < 2 is owing to the spin Casimir effect instead of the 
vanishing sublattice magnetization. Furthermore, our method provides an efficient and convenient 
tool that can estimate the correct exchange parameters and outline the quantum phase diagrams, 
which can be useful for experimental fitting processes in frustrated quantum magnets. 


I. INTRODUCTION 

Low-dimensional quantum antiferromagnetic systems 
have witnessed a great deal of interest for a long time 
due to their deep connection with magn etic properties 
of high-temperature superconductors .^321 However, low¬ 
dimensional quantum spin systems are of interest in their 
own right as fruitful resources of novel and exotic quan¬ 
tum phase^such as valence bond solidand spin liq¬ 
uids (SL).12EH11 The most indispensable ingredient in the 
emergence of these exotic quantum states is the quantum 
fluctuation caused by the SU (2) commutation relation 
of the spin operators. Frustration, on the other hand, 
acts as a very efficacious way of enhancing the quan¬ 
tum fluctuation effects and can even lead to the m eltin g 
of magnetic long range order at zero temperature.^^ If 
such magnetic long range order survives, it is expected 
that the quantum fluctuation effects will be suppressed 
and can have only small influences on the thermodynamic 
properties of the system.l^ Nevertheless, exceptions may 
occur, and as we shall see, quantum fluctuation can ex¬ 
hibit decisive consequences in non-collinear quantum an¬ 
tiferromagnets, despite the fact that the system is long 
range ordered. This can occur because of the Casimir 
effect generated by the zero-point fluctuation in a non- 
collinear background, which may generate some emergent 
phenomena to be less touched so far. 

The Casimir effect was originally discovered by Casimir 
in 1948, which states the presence of an attractive force 
bet ween two parallel conducting plates placed in the vac- 
uum.lSliSl This effect, which was described by Schwinger 
as one of the least intuitive consequences of quantum 
electrodynamics (QED), is actually ubiquitous in nature, 
covering many topics ranging from cosmology to con¬ 


densed matter physics jt arises when the quantum 
fluctuation of a general field (scalar, vector, spinor, or 
even tensor field) is modified by a static or slowly vary¬ 
ing ’’boundary”. This intriguing idea has generated con- 
tinui ng th eoretical interest in generalized Casimir prob¬ 
lems The same type of Casimir effects have been 
predicted and discussed in maw condensed matter sys¬ 
tems such as quantum liquid»ISl and nanoparticle sys¬ 
tems.^ The advantage of condensed matter systems as 
platforms to demonstrate the Casimir effects is the al¬ 
ready known structure of the quantum vacuum, at least 
in principle.^ Moreover, various exotic quantum phases 
in condensed matter systems may allow different char¬ 
acteristics of the Casimir effects. In this respect, the 
low-dimensional quantum magnets seem to provide an 
ideal playground for dealing with the generalized Casimir 
problem and a spin version of this intriguing effect may 
be expected consequently. 

The spin Casimir effect is the spin analog of the 
Casimir effect in vacuum, which describes various macro¬ 
scopic Casimir force and torque that emergent from 
quantum spin systems. Note the overall strength of the 
interaction generated by the Casimir Effect is propor¬ 
tional to the driving energy of quantum fluctuation (K) 
and its scale is related to the correlation strength of the 
fluctuations.lii^ Thus, the Casimir force is expected to be 
strong and long-ranged in a system with strong fluctua¬ 
tion and long range correlation. From this point of view, 
a system with highly degenerated ground states or in the 
vicinity of a quantum critical point is of particular in¬ 
terests and may generate rich Casimir physics. Such as 
the instability of charge ordered states caused by spin 
Casimir effect in doped antiferromagnets, where the zero 
point spin-wave fluctuation induce a uniformly attractive 
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force between hole clusters!^ Another interesting exam¬ 
ple is the quantum fluctuation lifted massive classical de¬ 
generacy of the ground st ate, w hich is called ’’quantum 
order by disorder” (QObD).E^ESlin this case, the ’’bound¬ 
ary” is the long range ordered classical spin structure. 
The effective description of the Casimir effect is appar¬ 
ent in some cases. For example, the Casimir (QObD) 
effect in a bilayer square-lattice Heisenberg antiferromag¬ 
netic model can be efficiently described by adding an ad¬ 
ditional term (S^ • to the original Hamiltonian.!^ 
However, there are other cases where the detailed form 
of the Casimir physics is obscure and can be seen only 
by loop expansions. 

In this article, we explore the emergence of spin 
Casimir effect in non-collinear ordered quantum mag¬ 
nets. The appearance of this effect is due to the zero- 
point fluctuation in a non-collinear ordered spin struc¬ 
ture and leads to the difference of measured ordering 
vector from the classical one. Although this difference 
has been discussed by several authors, its Casimir nature 
and related conse quences have not yet been thoroughly 
investigated! ^° ! ^^ ! This is mainly due to the fact that the 
spin Casimir effect is of order 0(1/5'), which makes the 
modification of the ordering vector much smaller than the 
classical value. Thus it is usually negligible in classical 
ordered systems. In contrast, we predict in this paper 
that in some circumstances a standard spin-wave the¬ 
ory becomes no longer applicable due to the presence of 
the spin Casimir effect, even though the system is long- 
range ordered. In this sense, the spin Casimir effect is 
no longer negligible. Furthermore, this effect can cause 
the spiral state instability, which is essentially different 
from other long-range order ’’melting” cases. We con¬ 
sider a two-dimensional spatially anisotropic triangular 
spin-1/2 antiferromagnet for the sake of general inter¬ 
est and perform a concrete and well-controlled calcula¬ 
tion. We believe that our results are equally applicable 
to other non-collinear ordered quantum systems with ar¬ 
bitrary spin value. 

It is known that an isotropic triangular lattice Heisen¬ 
berg antiferromagn et even for 5 = 1/2 may order into the 
so-called 120° state.!^* ^^^ * ^^ ! As the spatially anisotropic 
exchange interaction is turned on, the spin Casimir 
torque emerges, imposing modification to the classical 
ordering vector. Surprisingly, a careful 1/5 expansion 
in the anisotropic case shows that an usual perturbative 
estimation of the modification of the ordering vector be¬ 
comes divergent near the quantum critical point and the 
one-loop expansions of the energy spectrum and sublat¬ 
tice magnetization are strongly singular. These singular 
behaviors are believed to be the outcomes of the spin 
Casimir torque, and their appearance does not represent 
the onset of quantum disordered phases. To fix this point, 
we develop a self-consistent approach in the framework of 
the spin-wave theory, giving correct ordering vector mod¬ 
ification and excluding the singularities of the 1/5 expan¬ 
sions. Based on this self-consistent approach, a quantum 
phase diagram is obtained, which is qualitatively con- 
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FIG. 1: Exchange couplings between the different sites of the 
anisotropic triangular lattice. 


sistent with previous numerical works More than 

that, detailed results can be obtained with our approach 
by calculation that is no harder than a linear spin-wave 
expansion. Accordingly, our method can be considered 
as an efficient and fast experimental fitting tool for spiral 
phases. 

The remainder of this article is organized as follows. In 
Section H we introduce the anisotropic triangular lattice 
antiferromagnetic model, which is simple but instructive 
and grasps the core ingredients of physics. Section HI 
provides a brief review of the standard large-S expansion 
procedure and the formal definition of the spin Casimir 
torque which describes the modification of the classical 
spin structure due to quantum fluctuation. The first- 
order 0(1/5) quantum correction for the spin-wave spec¬ 
trum and its singular behavior are considered in Section 
IV A, while Section IV B is devoted to the calculation of 
the sublattice magnetization M to the order of 0(1/5^), 
which is divergent due to the presence of the spin Casimir 
torque. In Section V we develope the torque equilibrium 
spin wave theory (TESWT), which is free of the afore¬ 
mentioned singularities. Several physical properties are 
calculated within our approach and a quantum phase di¬ 
agram is obtained. And we discuss our scheme as an 
experimental exchange parameter fitting tool in Section 
VI. Finally, we draw our conclusions and discussions in 
Section VH. 


II. ANISOTROPIC TRIANGULAR LATTICE 
ANTIFERROMAGNET 


The two-dimensional triangular lattice is the simplest 
realization of geometrical frustration where a spin-liquid 
has been suggested. However, it is proved that the 
Heisenberg spins with the isotropic nearest-neighbor an¬ 
tiferromagnetic interacti ons on such lattice display an 
long-range ordered state.^^^^HlSl Nevertheless, the sub- 
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lattice magnetization is highly reduced from its classi¬ 
cal value due to the strong quantum fluctuation, indi¬ 
cating that small perturbations may destroy the long- 
range order an d driv e the system towards a quantum 
disordered stateP^^ In this respect, different kinds of 
interactions have been studied on the triangular lat¬ 
tice for the potential realization of the spin liquid 
state. Some of the most interested cases are the ring- 
exchange interaction, I 22 E 2 I the next-nearest-neighbor in- 
teracti orP ^^ ^ * ^^ !, and the spatial anisotropic interac- 
tion.IlIHlIIll] gg^gg jg particularly appreciated be¬ 

cause of its applicability to real mat erials such as inor¬ 
ganic and organic salts 

K - [BEDT - TT F) 2 Cu 2 iCNh and k - [BEDT - 
TTF)2Cu2[N{CN)\2^^ 

The Heisenberg antiferromagnet on an anisotropic tri¬ 
angular lattice has its Hamiltonian: 

<5i <52,(53 

?i = -Sj+J'^S,-Sj (1) 

(b> (b> 

where J is the interaction along the (5i, J' is the zig¬ 
zag interaction along the ^ 2,3 and the vectors 5i connect¬ 
ing neighboring sites are shown in Fig.l. In this work, 
both J and J' are positive and we denote the ratio of 
the coupling constants as a = J'/J. When J = J'^ 
the system is nothing but the isotropic triangular lattice 
Heisenberg antiferromagnetic model, whose ground state 
is a long-range ordered spiral state, i.e. the so called 
120° state. In the limit J = 0, the system is equivalent 
to the isotropic square lattice Heisenberg antiferromag¬ 
netic model, whose ground state is also long-range or¬ 
dered, the so-called Neel state. In the limit J' = 0, the 
system turns to the decoupled one-dimensional Heisen¬ 
berg antiferromagnetic chains, where long-range order is 
forbidden even at zero temperature due to the Mermin- 
Wagner-Coleman theorem.® In this case, the system is 
quantum disordered and shows many striking properties 
such as fractional excitation and power law correlation. 
Consequently, the related weakly coupled c hain region 
(J J') has attracted considerable attention.®2l23lMini 
We first sketch the classical case. The classical ground 
state is usually simple and constitutes the foundation of 
the further spin-wave expansion. In the classical case, 
the quantum fluctuation is absent and spins are vectors 
rather than operators. In the whole parameter space the 
ground state is a general spiral structure whose magne¬ 
tization Mi on lattice point is given by: 

riii = cos(Qc, • ri)x -I- sin(Q^, • ri)z (2) 

Here the spins are assumed to be in the x — z plane and 
the classical ordering vector = (Qci,0,0) is 

( 2 ' k , a > 2 

Qci = \ (3) 

t TT-b 2 arcsin(a/2), a<2 

which is determined by minimizing the classical ground 
state energy. This classical result is shown in Fig.2 and 
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FIG. 2: The ordering vector Q (in units of tt) of the opti¬ 
mal spiral state as a function of a = J'/J. The black curve 
is the LSWT result, the blue curve is the SE resultP^, the 
green curve is the MSWT resultp^, and the red curve is our 
torque equilibrium spin-wave theory (TESWT) result. The 
gray regions denote the values of a where the modified spin- 
wave calculations fails to converge indicating the onset of the 
quantum disordered phases.^ 



will be modified once the quantum fluctuation is consid¬ 
ered. 

Quite a number of theoretical approaches have been 
employed to treat the anisotropic triangular lattice 
Heisenbe rg m odel, such as linear spin -wave theory 
(LSWT),®® series expansions (SE),®® modified spin- 
wave theory (MSWT)I ® and density matrix renormal¬ 
ization group (DMRG).®® Besides, this model in the 
weakly coupled chain region has been studied by pertur¬ 
bative Bosoni zation and an effective Schrondinger equa¬ 
tion approach.®® The core prediction in the weakly cou¬ 
pled chain region is the existence of the long-sought-after 
spin liquid phases, which seems even more elusive after 
the variational quantum Monte Carlo (VMC) method 
predicted two spin liquid phases in this region.® Dif¬ 
ferent from the low a region, it is well recognized that 
the QObD effect stabilizes considerably the Neel state 
over the classical model, moving the Neel phase from the 
classical region a > 2 to a > 1.4, although the phase 
boundary determination is technique-dependent. How¬ 
ever, whether the quantum fluctuation spreads the tran¬ 
sition point between the Neel and spir al phases into a 
spin liquid is still controversial.®®*® At first glance, 
this controversy has nothing to do with the so called 
spin Casimir effect which only appears in the spiral phase 
around awl. On the contrary, the emergence of the spin 
Casimir torque in the spiral state naturally explains the 
instability of the spiral state, which provide the founda¬ 
tion of the QDbO effect and the potential existence of 
the quantum disorder phase, to be discussed below. 


































III. LARGE-S EXPANSION AND SPIN 
CASIMIR TORQUE 
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The spin-wave approach starts from a classical spin 
configuration which minimizes the Heisenberg interaction 
and treats the quantum de viati on from the ordered direc¬ 
tion as collection of bosons.^^^^In this work, the mapping 
from spin operators to bosons is performed via the Her- 
mitian Holstein-Primakoff transformation in a twisted 
frame.^ The resultant spin-wave Hamiltonian is given 
by 


'Htot — NEcI + lisw 



(a) 



(c) 



(b) 



(d) 


with 



=^2+^3 +^4 +0(5-1) (4) 

where Ed = 5^Jq is the classical ground state energy 
per spin, denote the terms of order 
extension is only to the cubic and quartic anharmonic 
terms. Such an approximation is sufficient for calculation 
of the 0(1/5) order result of spin-wave spectrum and the 
0(1/5^) order result of sublattice magnetization, which 
are our main interests. In the Fourier transformed repre¬ 
sentation, the explicit expression of various terms in the 
Hamiltonian reads as 

TD 

7/2 = 2S'^Aka|{^ak-^(akO-k + aj^alk) 

k 

773 = ( Ck(Qk-4-r)^^^P ^k^p^k+p) 

k,p 

7/4 = — 'y ( I [(Ai _3 -|- Ai -4 + A2-3 + A2-4) 

{k,} 

+ 77 i -4 + B2-3 -b B2-4) — {A4 + A2 

2 

-I-A3 + ^14)] 0^020304 • <^ 1 + 2 ,3+4 + 2 (-^l B2 

+B3){a[alala4 + 0402030 ^) • i5i+2+3.4| (5) 

Here, 1,2... denote ki,k 2 ..., and the following functions 
are introduced: 


Jk = J cos kx + 2 J' cos ^ cos -^ky 
2 2 

»7k = ^('^k-Q + dk+q), Ck = 2(5k-Q - Tk+q) 
^k = 2 ("^k + ??k — 2 Jq), Hk = -(Jk — ?7k) (6) 


On this basis, one may perform the 1/5 perturbation 
expansi on ei ther following the formalism developed by 
Belyaed2il25| o r turn ing to the Bogoliubovs quasi-particle 
representation.l221MI Here we follow the latter scheme, in 
which the quantum fluctuation induced singular behav¬ 
iors are much more evident. 

The quasi-particle representation is related to the 
Holstein-Primakoff representation by a Bogolyubov 


(e) (f) 

FIG. 3: The lowest-order cubic vertices [(a) and (b)], Hartree- 
Fock vertices [(c) and (d)] and ’’counterterms” vertices [(e) 
and (f)]. 


transformatiorPESE^ 

Ok = Uk&k + Wk&LkJ “k = “k&k + Wk&-k 


( 7 ) 


under conditions ^k ~ '^k = ^ 


2 I 2 ^k 

■“k + "Ck — ~ > 

Ek 


Bk 

2 ukVk = — 
^k 


( 8 ) 


with 


Ek = \ Al- Bl 


( 9 ) 


As a result, the linear spin-wave Hamiltonian takes the 
following diagonalized form: 

772=25^£k(&|,6k+J)-^ (10) 


and the cubic term 'H 3 in the new representation is 



( 11 ) 


The first term describes the magnon decay processes and 
is symmetric under permutation of two outgoing mo¬ 
menta. The second term serves as a magnon source and 
is symmetric under permutation of all three outgoing mo- 
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mentaP^ The explicit form of these two vertices are 


ri(l,2;3) — Ci^i(7273 + + C2^2(7i73 

+ K1K3) + C 3 K 3 ( 7 i 72 - K1K2) 


r2(l,2;3) — — Cl^l(7273 ~ «^2'«3) + C2K2(7 i73 


-K1K3) + C 3 K 3 ( 7 i 72 - K1K2) 


( 12 ) 


with 


? = \/£i£ 2£3, Hi = y/Ai + Bi, 7i = ^JAi- Bi (13) 
where i G (1, 2, 3). 

It is noted that the transformation of the quartic terms 
is rather cumbersome. Given that we are only inter¬ 
ested in the one-loop results, the quartic terms can be 
conveni ently decoupled using the Hartree-Fock approx- 
imation.ElB3H26| Introducing the following Hartree-Fock 
averages in momentum space 


Wk = (aiLak) = ^ Ak = (aka_k) = ^ (14) 


A 


1 


Bk 


2 £k 


2 £k 


the quartic terms turn to the form 7^4 = Ei + 6 'H 2 , where 
E 4 is the Hartee-Fock corrections to the ground state 
energy and 6 'H 2 is the 1/S modification to the harmonic 
spin-wave Hamiltonian with the form 


Sn2 = Y. Sskblbk - ^{bkb^k + blbl^) 


(15) 


where 


fek = (wk + i^k)<5^k - 2ukVkSBk 
Ok = (wk + vl)SBk - 2ukVkSAk 


(16) 


and 

JHk = Ak 


1 r 


7lp(7lk—p — Ak — Ap — 77k—p) 


E 

p 

+Bp{-^ + Bp ^ 

SBk = 77k — ^— 77p(Hk-p-^ — Hp — 77k-p) 

p ^ 


+^p(77k + Bp) 


(17) 


The effective Hamiltonian that combines all the terms 
together now reads 

=Y [( 2 * 5 ek + 5£k)bibk - ^{bkb-k + blbl^) 


k 

+i\ 

1 




k.p 


+ ^r2(p,-k-p;k)6j,6l^_p5T -H.c. 


(18) 


which provides the basis for the 1/S perturbative ex¬ 
pansion. The related Feynman diagrams are shown in 
Fig.3(a)-(d). 

At the same time, the 1/S expansion contributes to the 
corrections of the ground state energy as well. This modi¬ 
fication comes from the zero-point fluctuation of magnon, 
which is the fluctuating vacuum of our Casimir problem. 
For the sake of simplicity, we only consider the first order 
corrections to the vacuum energy per spin 

Eyac = E,i +E 2 = S^Jq + 5( Jq + ^ £k) (19) 

k 

Here E 2 is the energy correction from the harmonic spin- 
wave fluctuation. 

With this vacuum energy modification, the ordering 
vector of the system should be determined by minimiz¬ 
ing the modified vacuum energy Eyac via SEyac/SQ = 0. 
However, this simple variational equation can not be 
solved directly due to the fact that the spin-wave spec¬ 
trum function is only well-defined at Q = Q^,;. Thus, 
the variation is normally treated approximately as an ex¬ 
pansion around Q^;. The 1/S order result is 

Q = Qci + Qi (20) 


with 

Qi 


1 


Ak + 77k 

9Jk+Q 

2S 


A 

dQ 


( 21 ) 


These ordering vectors have only the x component in our 
case via = ((3i,0,0) and Q = ((5,0,0). This result 
seems reasonable and is usua lly treated as the new order¬ 
ing vector of the system.l2212U However, this is not the case 
when the system is in the vicinity of a quantum critical 
point. The first order correction Qi goes to infinity as a 
approaches the classical spiral/Neel critical point a = 2, 
which is shown in Fig.4. 

Clearly, this result is unphysical and needs to be regu¬ 
larized. Here we propose a torque description of the spin 
Casimir effect, which is analogous with the force descrip¬ 
tion of the conventional Casimir problems in QED .1141121 
To access a quantitative description, we define the spin 
Casimir torque as 


T«c(Q)=^ 

k 



Ik yac 


( 22 ) 


where iT^ac) represents the quasi-particle vacuum state. 
Notice that Esc is a function of Q defined on bonds and 
represents the tendency of the modification to the rel¬ 
ative orientation of each spin. One can have the same 
definition for a QObD system, which is equivalent with 
other conventional description methods. The main ad¬ 
vantage of the torque description is to provide a much 
more intuitive picture. As an example, the resultant 1 /S 
order spin Casimir torque at is 


^ _ -S' Ak -b 77k dJk+Q 


Qc 


(23) 
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FIG. 4: Tsc, Qi/tt and F(Q-) as a function of a. Here Q- = 
Q — 0 ^. The insert is these functions in region 0 < a < 1. 
Here both Qxjix and F{Q-) are divergent when approaching 
a = 2 as described in the main text. 


And in our case Tsc(Qci) has only one component noted 
as TsciQci)- Clearly, this spin Casimir torque is well- 
defined in the whole parameter space as shown in Fig.4. 
And the zero ordering vector modification at a = 1 is eas¬ 
ily explained by the zero spin Casimir torque due to the 
triangular symmetry. For a ^ 1, the spin Casimir torque 
tends to arrange the spins connected by the strongest 
bonds collinearly. This tendency is consistent with the 
conventional statement that fluctuation favors collinear 
spins. Based on the torque interpretation and noting 
that S^d'^jQ/dQ'^ is just the classical spin stiffness ps, 
the complex Eq.(21) is nothing but the Hooke’s law for 
spin system 

TsciQci) = —PsQi (24) 

Additionally, the reason of the divergence of Qi in ap¬ 
proaching the critical point {a = 2) is also apparent. The 
transition at a = 2 is continuous and the spin-wave ve¬ 
locity for Goldstone excitation vanishes at the transition. 
When the system approaches a = 2, the quantum fluc- 
tuat ion d ominates and the spin stiffness quickly goes to 
zeropSMI Although Tgc also goes to zero as a —f 2, its 
Casimir nature is enhanced as fluctuation dominates. As 
a consequence, it approaches zero much more slowly than 
Ps, and eventually the resultant Qi blows up around the 
critical point, exhibiting the singularity at the point. One 
may suggest that this singular behavior of Qi is under¬ 
standable because the divergence of the fluctuations near 


a Lifshitz point is well known.l^ As a consequence, the 
spin-wave th eory should be invalid when near a quantum 
critical point .1^ However, the difficulty for the spin-wave 
theory is the existence of other singularities even far away 
from the critical point, in the presence of the spin Casimir 
torque. 


IV. SPIN CASIMIR EFFECT INDUCED 
SPIN-WAVE SINGULARITY 

The core ingredient of the spin-wave theory is the ex¬ 
pansion in powers of 1/S around the classical saddle 
point. Strictly, this the ory is only well-defined in the 
large-S limit {S ^ and thus less effective for low¬ 
dimensional quantum spin systems, as quantum spin fluc¬ 
tuations typically increase in reduced space dimensions 
and for small spin quantum numbers S. It is surprising 
to observe that the standard spin-wave approach can give 
very accurate description of the zero-temperature physics 
of a number of low-dimensional spin models such as the 
S' = 1/2 Heisen berg ant iferromagnets on square and tri¬ 
angular lattices.III^I 2 iMl expansion ap¬ 

proach can still be considered as a useful technique if an 
ordered state is observed. Nevertheless, it fails once the 
spin Casimir effect is taken into account, to be demon¬ 
strated here. 

In the subsequent two subsections we will show that 
the one-loop expansions of the energy spectrum and sub¬ 
lattice magnetization are strongly singular and these sin¬ 
gular behaviors are related to the spin Casimir torque. To 
access the explicit structure of these singular behaviors 
and show the breakdown of the conventional 1/S expan¬ 
sion procedure, we first ignore the spin Casimir effect and 
Q is identified as Q^;. 


A. Spin-wave Spectrum 

Perturbative expansion for the spin-wave spectrum has 
to be performed order by order in 1/S and takes into ac- 
count all the quantum corrections of the same order.l^^ 
Only in this manner one can ensure cancellation of all 
possible divergences in the individual contributions and 
preserve the Goldstone theorem.l^ The first order cor¬ 
rection is straightforward. The new pole of the magnon 
Green’s function is determined by the so called Dyson 
equatioiP^ 


e = Ek ■ 


-I-— + S3(k, e) -I- S3(k,e) (25) 


with 


S^(k,e) 


1 |ri(p; k)p 

2 ^ - £p - £k-p + i0+ 

p ^ ^ 

_ 1 |r2(p; k)p 

2 £ -|- £p -|- £k+p — *0"*" 


( 26 ) 
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FIG. 5: The lowest-order normal [(a) and (b)] and anomalous 
[(c) and (d)] magnon self-energies generated by cubic terms. 


The diagrammatic representations of the normal self¬ 
energies from the cubic terms are shown in Fig. 5 (a) and 
5(b). 

This equation can be solved either self-consistently 
through the off-shell approximation or by replacing e 
with linear spin-wave spectrum ek, i.e. the so called 
on-sell approximation. The 1/5 order correction to the 
spectrum F'(k) is obtained within the on-shell approxi¬ 
mation, which leads to the following expression for the 
renormalized spectrum: 


F"(k) — Ek + 


1 


fek + S3(k,£k) + S3(k,£k) 


(27) 


with 


Go(k) 

GQ(k) 

Lo(k) 

LqO^) 


P 

V- 7p7Lp 


___ t^2 


Sn r 


E 


r2 


£p + ^k- 


(29) 


and 


ro = 

rQ = 


1 

1 


Ck(^p^k—p 


[Cp^p7k—p 


7p7k—p) (Cp “b Ck—p)^p^k—p 
Ck—p7p^k—p] 


These equations allow a straightforward examination of 
the Goldstone mode at k = 0 and k = Q. Noting that 
7 o = 0 and Co = 0, it is easy to prove Go(0)/eo = 0 and 
Lo(0) = 0. Hence, as fc —>■ 0 we have 


F(0) 


25 


Go(k) 

2ek 


^o(k) 

4 


k^O 


= 0 


(31) 


Indeed, the Goldstone mode is preserved at k = 0. 

In contrast to the k = 0 case, the Goldstone mode 
usually appears by cancelation among several terms.123^ 
Notice that kq = 0 we obtain 


F(Q) 


A 

25 


gQ(k) 

2ek 


Lqjk) 

4 


k^-Q 


(32) 


Based on this expression, the one-loop spin-wave spec¬ 
trum F"(k) can be easily obtained by numerical integra¬ 
tion of self-energies. However, the numerical result shows 
the singularity of the spectrum at k = Q (as shown in 
Fig.6.) and absence of the Goldstone mode, while the 
Goldstone mode is usu ally expec ted at every order of the 
perturbative expansionl ^ l ^'^ l ^^ l S H Thus the absence of the 
Goldstone mode and the instead appearance of the sin¬ 
gular behavior are quite unusual. 

To understand this singularity, a careful examination 
on all the contributions is needed. Regarding the Gold¬ 
stone excitations, the terms proportional to Ek can be 
ignored safely. The resultant explicit form of the self¬ 
energies in the expression of F'(k) is 


In this case, Gq(Q) and Hq{Q) do not equal zero in an 
obvious way, the resultant self-energies are divergent at 
k = Q. To get the exact analytic structure, we introduce 
the following useful equalities 

Ck — '^k-Q ~ 7k 
Cp = Kk-p “ 7p + ^1 
Ck-p = Kp“7Lp+^2 

I 2 2 

<^1 — '^p-Q ~ ^k-p 

<^2 = '^k-p-Q ~ ^p (33) 

noting that i5i and <52 equal zero at k = Q. The ex¬ 
pansion of F{Q) is straightforward with the aid of these 
equalities, and the final result is 


(5Ek 




^kG'o(k) + 7kGQ(k) 
4^o(k) -b7k-^Q(k) 


( 28 ) 


F{Q) 


Aq Tip + Bp ^Jp-^-Q 

Z5 ^ ^ dQ 


k- Q 


Ek 


k=Q 


^T,,(Q) . 

S (Sk 


k=Q 


( 34 ) 








































FIG. 6: Numerical results of the spin-wave spectrum around 
k — Q for a = 0.5. 


The explicit expression along fcy = 0 is 


F[Q) = 

VQ „ £p 


= sgn(fca; - Q)^—Tsc(Q) 

J VQ 


(35) 


Here sgn(a:) represents the sign function and vq is the 
spin-wave velocity = Q along ky = 0. The an¬ 

alytic result shows that our spin-wave spectrum behave 
like a step across k = Q rather than a cone, which is con¬ 
sistent with our numerical result. More seriously, F(Q) 
also blows up as a approaches the classical critical point 
q; = 2 as shown in Fig.4. 

Interestingly, this singular behavior is directly con¬ 
nected with the spin Casimir torque defined previously. 
It is indeed due to the ignorance of the spin Casimir effect 
in the above treatment. Since the spin Casimir torque is 
directly connected with the 1/S ordering vector correc¬ 
tion Qi, one suggests that the Qi induced modification 
of the linear spin-wave spectrum Sk may avoid the above 
problem in some region of the parameter space where Qi 
is still finite and reasonable. In fact, with the equality 

VQ = 2^PsAq (36) 


the relationship between F{Q) and Qi is apparently 
shown as: 


sgn{k^ -Q) 

P[Q) = -- VQ ■ Qi (37) 

On the other hand, a direct I /S order correction due to 
the ordering vector modification is 


SF{Q) 


d£k 

dQ 


■ Qi 


(38) 


which does not match our expectation. In other words, 
this singular behavior can’t be regularized in an conven¬ 
tional 1/S manner, and it can even lead to divergence in 
high order expansions. 


B. Sublattice Magnetization 

In this section, we turn to the sublattice magnetiza¬ 
tion which is the order parameter for general long-range 
ordered spin states. It can be used to probe the possi¬ 
ble existence of quantum disordered phases. Within the 
spin-wave approach, its definition is 

(S) = S - (aja,) = S - SSi - SS 2 (39) 

Here the first quantum correction SSi is the linear spin- 
wave result given by 


SSi=J2nk (40) 

k 

And the second term is the second-order correction to the 
sublattice magnetization, which needs an evaluation of 
the one-l oop re sults of the normal and anomalous Green’s 
functions This term have three contributions 


SS 2 = SSi + SSi + SSi (41) 


These contributions are related to the normal and 
anomalous self-energies as shown in Fig.3(c), 3(d) and 
Fig.5. The contributions from the normal self-energies 
are collected as SSj, and those from the anomalous self¬ 
energies Ok and Eg'’* are collected as SS 2 , 6 S 2 respec¬ 
tively. Calculations of each contribution to the second- 
order correction are straightforward and the final results 
are 


^^2 = ^ E (J") ’ ^^2 = ^ E (>") 


k k 

with i e (1, 2, 3), ItotiM) = and 

|r2(k,p)p 


^i(k) 

^2(k) 

^3(k) 


1 Ak 

9 TT ' 


2 ^k (^k T T ^k-l-p)^ 


1 Hk Ok 

1 . ri(k,p)r 2 (-k,p) 

2 p ^k(^k 4“ ^p 4“ ^k—p) 


(42) 


(43) 


As often taken for the high-order perturbative expan¬ 
sion, the integrands of all the three contributions are di¬ 
vergent. Especially, the integrands / 2 (k) and / 3 (k) be¬ 
have as 0(1/A:^) at k —>■ Q, implying that not only the 
leading divergences in them, but also the subleading ones 
0(I/fc^) must cancel in order to produce finite result. 
Expanding near k = 0 and k = Q points, such a can¬ 
celation can be verified analytically at a = 1 as shown 


k=Q 
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by Chubukov!^ Here we perforin the expansion with a 
general a. 

Because the divergence appears only near k = 0 and 
k = Q, the terms proportional to Ek (finite part) can 
be ignored as before. After lengthy but straightforward 
derivation one obtains 


Green’s function, but it is generated from the abnormal 
Green’s function rather than the normal Green’s function 
which accounts for the spin-wave spectrum. 

Correspondingly, a conventional 1/S consideration of 
the ordering vector modification induced sublattice mag¬ 
netization correction is^ 


with 


/i(k) 

/2(k) 

/3(k) 


8Ek 

4e3 

8 el 


K^ 17 o(k) -I- 7 ^t 7 Q(k) 
-^2Go(k)+72GQ(k) 
KlVo{k)-^lVQ{k) 


Uo{k) 

Uqik) 

K)(k) 

Vqik) 


E 

P 

E 

p 

E 

p 

E 


r2 
^ 0 


(sk + + £k-p)^ 

_ 

(^k ^p ^k—p) 

_)_o_ 

£k + £p + ^k-p 
Ek ~t“ Ep -f Ek—p 


(44) 


(45) 


where Go(k), GQ(k), Fq, and Fg have been defined in 
Eq.(29) and Eq.(30). These equations are similar to those 
obtained in the expansion of spin-wave spectrum. As a 
consequence, the divergence cancelation results are quite 
similar as well. The leading and subleading parts of the 
integrands are zero at k = 0 point, as expected. How¬ 
ever, the perfect divergent cancelation near k = Q is 
not accessed in our general case and the final result is a 
subleading divergent contribution 


/*"(Q) = ^p + ^p ^^p+Q k-Q 


2e2 

AqHq 


T5c(Q) 


Jp dQ 

k-Q 


Ek 


Ek 


k=Q 

(46) 


k=Q 


whose explicit expression along fcy = 0 is 


iTotiQ) = sgn(fc,-Q) 


Hq SAq 

VQ 


E 


Ap + Hp dJ, 


p+Q 


= sgn(fca; - • Tsc(Q) 




= ■ F{Q) 


dQ 


(47) 


which is related to the spin Caisimir torque too. One 
may argue that this divergent result is inherited from the 
singular behavior of the spin-wave spectrum. As a matter 
of fact, this divergence is indeed caused by the one-loop 


< = (48) 

^ Qci 

which is a finite integration. The total sublattice mag¬ 
netization is still divergent even by adding this contribu¬ 
tion, preventing the spin-wave expansion beyond the har¬ 
monic approximation in spiral phases. As we have shown, 
the only region in the parameter space that is free of di¬ 
vergence is the a = I point, but it is not likely that the 
long-range order only exists at this single point .l23EIlpui-_ 
thermore, one can expect more serious divergence when 
the expansion to higher order is carried out, making the 
conventional spin-wave theory failed. Accordingly, an al¬ 
ternative expansion scheme that can encompass the spin 
Gasimir effect is needed. 


V. TORQUE EQUILIBRIUM SPIN WAVE 
THEORY 


The spin Gasimir effect described in our work can be 
generalized as the effect to shift the quantum fluctua¬ 
tion induced classical saddle point in quantum spin sys¬ 
tems. As a consequence, the corresponding spin exci¬ 
tations should be considered by expansion around the 
shifted (new) saddle point. However, a shift of the sad¬ 
dle point is prevented by two essential issues: 1. Where is 
the new saddle point while the 1/S expansion results are 
divergent at some parameter region? 2. How to shift the 
saddle point from to Q in the effective Hamiltonian 
with so many functions defined only at Qj,;? 

The first issue can be easily handled using our torque 
description of the spin Gasimir effect. The quantum fluc¬ 
tuation induced shift of the saddle point is accomplished 
by the spin Gasimir torque. As this spin Gasimir torque 
shifts the spin structure away from its classical saddle 
point, a classical reaction spin torque is generated by the 
deformation. The final saddle point is determined by the 
torque equilibrium condition 


T,e(Q)+T,((Q) =0 


where 


Tei(Q) 


dE,i{Q) ^ 

5Q aq 


(49) 


(50) 


This is nothing but the variational equation of the vac¬ 
uum energy SEyac/^Q — 0. However, this equation can 
not be solved directly owing to the existence of the sec¬ 
ond issue, which is the main difficulty in our perturbative 
expansion scheme. To overcome it, an alternative spin- 
wave expansion scheme is needed. 
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It is noted that and i3k are well-defined for all 
Q. The ordering vector can be shifted arbitrarily for the 
Hamiltonian 7 ^2, 7^3 and in Eq.(5). Yet, for Q ^ 
the resultant ek(Q) would be imaginary at some k points 
reminding that the expansion is carried out around the 
wrong saddle point. The basic idea of our scheme is to re¬ 
arrange the perturbation processes and take into account 
the ordering vector modification in an self-consistent way. 

Given Q as the final ordering vector of the system, all 
the functions in Eq.(5) are shifted from to Q. Eor 
any function, taking Hk as an example here, it can be 
written as 

^k(Q) = Ik(Q) + A£(Q) (51) 



0.0 ^ I W -- atJ 

(0,0) (rO) (2r0) 


with 


H£(Q) = Hk(Q)-Hk(Q) (52) 

where Hk is the Hk function of another spin system whose 
classical ordering vector equals to Q. Here we assume 
that this spin system has the same symmetry and set 
of exchange integrals as the original one. As a result, 
this spin system is nothing but the original one with a 
different a denoted as 5. The spin-wave Hamiltonian 
written in this manner is 

•H2(a,Q) = H2(5,Q)-fH^ (53) 


with 

^ _ _ TD 

7^2 = 25^Akaj^ak-^(ak^-k + aiLaLk) 

k 

'^2 = ^^^Akaj^ttk-^(aka-k + aj^a^k) 

k 

The 772 is free of imaginary energy problem and 77^ is 
obviously proportional to In spite of the 25 factor, 
the Hamiltonian 77^ is of order 0(5°) just as 774 due to 
the fact that A^(Q) and 77^(Q) are of order 0(1/5). In 
the spirit of the 1/5 expansion, we treat 77| as an interac¬ 
tion term and 772 as the modified harmonic Hamiltonian. 
The spin Casimir torque at the equilibrium point can be 
expressed approximately as 

k 

This approximate spin Casimir torque has only one com¬ 
ponent TgciQ) as well. And the torque equilibrium equa¬ 
tion can be approximately written as 

QJq _ 1 Ak + Hk cIJk+Q 

k 

The exchange parameters on the left hand side of the 
equation is exact as a while the parameters on the 
right hand side approximate as a. Noticing that a = 


FIG. 7: Spin-wave spectrum for a — 0.5 along the symmetric 
directions in the Brillouin zone. The blue line is the results 
obtained in the LSWT approximation. The red and green 
lines are the real (energy) and imaginary (damping rate) parts 
of our 1/5 results, respectively. And the gray areas show the 
width of the spectral peaks due to the damping.^ 


—2cos(Q/2), this equation can be solved and the numer¬ 
ical results are shown in Fig. (2). The results show drastic 
modification of the classical ordering vector caused by the 
quantum fluctuation. And our results are similar with 
the results o btained by much more sophisticated numeri¬ 
cal methods.l^^ISIESI Interestingly, for a > 1.2, the torque 
equilibrium equation has no non-trivial solutions other 
that Q = 2tt, hence no spiral order is stable. Therefore, 
in this region, the system is either in a QObD induced 
Neel phase or some quantum disordered phases. Never¬ 
theless, whether a spin liq uid ph ase exists in this region 
is still under controversial.l2Zl22JI32l Either way, the con¬ 
ventional spin-wave expansion starting from the classical 
spiral state is qualitatively incorrect for 1.2 < a < 2. As 
5 becomes large, the range of this spiral instable region 
is narrowed and eventually shrinks to a point a = 2 for 
5 = oo. 

The above separation procedure can be repeated to ev¬ 
ery term in Eq. (5). Eventually, we obtain a whole set of 
77i and 77/. The resultant new terms are very similar 
to the counterterms in quantum field theory.l^ Interest¬ 
ingly, the Casimir problem in the quantum field theory 
is divergent and this divergence can be regulated by in¬ 
troducing the counter-term s tha t are fixed with physical 
renormalization conditionsIt is expected that our 
divergence problem can be solved in a similar way. Given 
that we are only interested in the results at 1/5 order, 
the ’’counter-terms” 77§ and 77/ can be neglected. 


77«„, = 772 + 77/+ 773 + 774 (57) 


Following the same procedure described in Section HI, 
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the effective Hamiltonian reads 


Ueff = { (25?k + - ^(&k6-k + 


+25 


4^k^k - ^{bkb-k + 


+iV^Yl [^ri(p,k-p;k)MLp^P 

k.p 

+ ^r 2 (p,-k-p;k)&j 36 Lk-p^k-H.c. 

Here R represents i?(5, Q) and 


(58) 


4 = (Wk + ^)^k - 2uk^'k-Bk 
^k “ ('*^k "f '^k)'®k ~ 2MkWk^k (59) 


The diagram representation of these extra terms are sim¬ 
ilar to the counter-terms in the quantum field theory as 
shown in Fig.3(e) and 3(f). The resultant contribution 
from our ” counter-terms” to the 1/5 spin-wave spectrum 
is 


^k = 


1 - 
Ek - 


^k^k(Q) — ^kSk(Q) 


- Ek 


(60) 


whose contribution to k = 0 point is zero while that to 
k = Q point is 




dJQ k- Q 
dQ Ek 


(61) 


k=Q 


And the contribution from the other terms is 


F(Q) = ^T,,(Q) • 

5^ Ek 


(62) 


k=Q 


Noting that the approximate torque equilibrium equation 
Eq.(56) can be written as 5 ^i9Jq/ 9Q = —Tsc(Q), the 
singular contribution is canceled by the ” counter-terms” 
contribution. The cancelation along other directions is 
straightforward to certificate and the Goldstone theorem 
is preserved as it should be. The numerical results of 
the spin-wave spectrum are shown in Fig.7, which man¬ 
ifests well b ehaved Goldstone modes and magnon decay 
effects.l^^lSlHISIxiie downward modification of the magnon 
spectrum and the line shape of the magn on damping rate 
are consistent with previous work^^filMl, Por a = 1, the 
spin Gasimir effect vanishes and our results are exactly 
the sa me with that obtained by Ghernyshev and Zhito- 

mirsky.l2SlM155l 

On the other hand, the contribution to the sublattice 
magnetization from our ’’counter-terms” is 

= ( 63 ) 



FIG. 8: The sublattice magnetization (5) as a function of a = 
J'/J. The red triangular dots represent the LSWT results, 
and the purple red dots represent the standard 1/S expansion 
results. The blue square dots are our torque equilibrium linear 
spin-wave theory (TELSWT) results, and the green red dots 
are our 1/S-TESWT expansion results. 


with 

^c(k) = ^|^-^ (64) 

^ Ek Ej^ 

This integrand is also zero at k = 0 point and divergent 
at k = Q point as 




Q 


dQ 


Ek 


(65) 


k=Q 


And the divergent contribution from the other terms are 




^qRq 


T5c(Q) 


k- Q 


Ek 


k=Q 


( 66 ) 


These two divergences cancel again owing to the approx¬ 
imate torque equilibrium equation and this cancelation 
persists to other directions. Eventually, we can obtain a 
finite second order correction to the sublattice magneti¬ 
zation. 

The above results allow us to perform systematic cal¬ 
culations on the sublattice magnetization of the system 
in the whole parameter space, which may indicate po¬ 
tential existence of spin liquid phases.l^ The calculation 
is straightforward and follows our previous established 
scheme in the spiral region and the sublattice magne¬ 
tization in the Neel state is obtained within a standard 
spin-wave expansion scheme. The final results of the sub¬ 
lattice magnetization and related phase diagram are pre¬ 
sented in Fig.8. In the obtained quantum phase diagram, 
the spiral state is destroyed by quantum fluctuation as 
a < 0.5, which is not far from the MSW results. In addi¬ 
tion, above the isotropic point a = 1, the spiral state is no 
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longer stable in the region a > 1.2, which is also consis¬ 
tent with previous numerical studies.^^lllllxhis instability 
is not due to the vanishing sublattice magnetization but 
owing to the fact that the classical saddle surface is so 
shallow that quantum fluctuation can induce the instabil¬ 
ity of the classical spiral state. In our torque description, 
this instability can be understood as the so weak stiffness 
of the ’’spring” that the spin Casimir torque can squeeze 
it arbitrarily until the torque disappears. In this man¬ 
ner, the spin Casimir effect shows us a spatial way to 
’’melt” a long-range non-collinear ordered state, which is 
different with the perception of the collinear cases. What 
should be mentioned here is that the exact structure of 
the phase di agram in the region 1.2 < a < 2 is still un- 
ciear ,121130131] and the perturbative nature of our approach 
prevents us to identify whether the system is in a QObD 
induced Neel phase or some quantum disordered phases. 

Our approach can be easily generalized to systems with 
multi parameters by considering more counter-terms such 
as ?{§ and Other than that, more self-consistent 

equations can be obtained using the divergence cancela¬ 
tion condition in higher order. The non-renormalizable 
form of the effective Hamiltonian ensures that we can 
have all kinds of divergences and counter-terms as long 
as we carry out our calculations to high enough order.l^ 
However, the whole spin-wave basis can breakdown if 
high order loops are considered and it is usually unnec- 
essary.E One can simply apply the present scheme to the 
parameter dominating the ordering vector, given other 
parameters unchanged. In addition, our approach can 
be considered as a general approximation method for the 
quantum fluctuation induced saddle point shift problems. 


VI. EXPERIMENTAL APPLICATIONS 

Due to the fundamental theoretical interest, further 
impetus to investigate frustrated magnets has arisen from 
recent experimental developments identifying several ma¬ 
terials in w hich s pin liquid like unconventional behaviors 
is observed Hence potential progress is expected in 
the near future. So as to test the theory experimentally, 
a precise information on the spin-Hamiltonian parame¬ 
ters for the materials of interest is demanded. And the 
most effective way of solving this problem is to suppress 
quantum fluctuation by strong-enough magnetic field. 
However, such experiments are only possible for those 
systems with small-enough exchange parameters such as 
CS 2 CUCI 4 , so that the required field can be available 
practically. It would be highly desirable to dispose of a 
fast tool that can estimate the correct exchange param¬ 
eters and outline the quantum phase diagrams. 

Here we propose that the methods presented in this 
work will serve this very purpose. Compared with other 
sophisticated numerical and analytical methods, our ap¬ 
proach is much less technique involved. As a matter of 
fact, detailed results can be obtained in our approach 
by calculation that is no harder than a linear spin-wave 


expansion. Other than that, our approach can be con¬ 
veniently applied to magnetic systems with general spin 
Hamiltonian and the only requirement is the existence of 
the long-range order. In other words, our scheme only re¬ 
quires experimental results with magnetic fields that can 
drive the system into long-range ordered states. More 
importantly, the long-range ordered phases in magnetic 
fields are usually non-collinear, such as ’’fan” and ’’um¬ 
brella” phases. Add it all up, our method shows clinical 
improvement compared with the conventional parameter 
fitting processes based on LSWT but keeps the compu¬ 
tational difficulty nearly unchanged. 

In order to demonstrate the superiority of our method, 
we consider CS 2 CUCI 4 as an example. The effective ex¬ 
change parameters obtained from the global fit of the 
zero-field neutron scattering results are 0 . = 0.175.1^ This 
fitting result based on LSWT is far from the exact result 
a = 0.34 that obtained using the high-field tech nique 
where quantum fluctuations are quenched out How¬ 
ever, by using our approximate torque equilibrium equa¬ 
tion Eq.(56), the LSWT fitting result 5 = 0.175 predicts 
the bare parameters with a = 0.35, which is very close 
to the high-field exact result. Beyond that, our approach 
further shows that the sublattice magnetization of the 
system without Dzyaloshinskii-Moriya (DM) interaction 
vanishes based on linear approximation. As a result, the 
long-range spiral state in CS 2 CUCI 4 can’t be stable with¬ 
out DM intera ction, which is consistent with other theo¬ 
retical analysisl ^^ l ^^ l 

In summary, our approach enables one to estimate the 
correct exchange parameters and outline the quantum 
phase diagrams based on simple linear spin-wave approxi¬ 
mation. Thus, besides the theoretical setups, our method 
can also serve as an efficient tool for the experimental 
fitting processes of the exchange parameters in general 
frustrated quantum magnets. 


VII. SUMMARY 

We have presented a detailed analysis of the spin-wave 
expansion on the spatial a nisotro pic triangular lattice 

Heisenberg antiferromagnets.^^HlII phenomena that 

the classical ordering vector is modified by the quantum 
fluctuation is carefully studied and its Casimir nature is 
revealed. This effect sh ares t he same origin with the 
well-known QObD effect Both of these cases can 
be interpreted as the spin Casimir effects in which the 
quasi-particle vacuum with zero-point fluctuation plays 
the role of the fluctuating vacuum and the classical spin 
structure act as the macroscopic boundary. To describe 
these spin Casimir effects quantitatively we further de¬ 
fine a spin Casimir torque, which describes a long-range 
torque effect generated by quantum spin fluctuation. 

Base on these results, it is shown that the presence of 
the spin Casimir effect can induce divergent results in a 
conventional spin-wave expansion even though the long- 
range order is stable. A careful expansion shows these 
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divergences are directly connected with the spin Casimir 
torque. As a result, the appearance of these divergences 
inv alida te the conventional 1/S expansion in an obvious 
way.li^ And the encountered problems are rather generic 
and common to a variety of frustrated antiferromagnets 
regardless of the spin value. For the systems with large 
spins, the situation is especially aggravating even though 
the spin Casimir torque goes to zero as S' —> oo. This is 
due to the fact that long-range order is more stable and 
the 1/S expansion turns to be more reliable when S 3> 1 
while the divergence prevent any reasonable prediction. 

In the present work, we have developed a self- 
consistent approach in the frame of the spin-wave the¬ 
ory, which is applicable to variety of systems with quan¬ 
tum fluctuation induced saddle point shift problems. The 
self-consistently calculated modification to the ordering 
vector is finite and close to previous SE and MSW re¬ 
sults as shown in Fig.(2). Furthermore, our approach 
regularizes all the divergences in the 1/S expansion ef¬ 
fectively. This accomplishment allows us to calculate 
many spin-wave properties beyond linear approximation. 
The spin-wave spectrum results are shown in Fig. 7, which 


present w ell beh aved Goldstone modes and mangnon de¬ 
cay effects.1^^^2®^ Other than the conventional spin-wave 
properties, an approximate quantum phase diagram is 
also obtained, which d isplay s good consistency with pre¬ 
vious numerical works.^^^^ These results evidence that 
our approach is a suitable tool to study various other 
problem in non-collinear quantum antiferromagnets. Be¬ 
sides the theoretical setups, our method can be further 
applied to estimate the correct exchange parameters and 
outline the quantum phase diagrams, which can be useful 
for experimental fitting processes in frustrated quantum 
magnets. 
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